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We present the first direct measurement of the gravity-field curvature based on three conjugated 
atom interferometers. Three atomic clouds launched in the vertical direction are simultaneously 
interrogated by the same atom interferometry sequence and used to probe the gravity held at three 
equally spaced positions. The vertical component of the gravity-held curvature generated by nearby 
source masses is measured from the difference between adjacent gravity gradient values. Curvature 
measurements are of interest in geodesy studies and for the validation of gravitational models of the 
surrounding environment. The possibility of using such a scheme for a new determination of the 
Newtonian constant of gravity is also discussed. 


In the last two decades, atom interferometry jl] has 
profoundly changed precision inertial sensing, leading 
to major advances in metrology and fundamental and 
applied physics. The outstanding stability and accu¬ 
racy levels SI combined with the possibility of easily 
implementing new measurement schemes |^-[7| are the 
main reasons for the rapid progress of these instruments. 
Matter-wave interferometry has been successfully used 
to measure local gravity Q, gravity gradient I9H11|. the 
Sag nac effect JT^l , the Newtonian gravitational constant 
fl3l - {T6l ]. the fine structure constant 0 , and for tests of 
general relativity 00. Accelerometers based on atom 
interferometry have been developed for many practical 
applications including geodesy, geophysics, engineering 
prospecting, and inertial navigation 20 h 22|. Instruments 
for space-based research are being conceived for differ¬ 
ent applications ranging from weak equivalence principle 
tests and gravitational-wave detection to geodesy [23I . 0 . 

One of the most attractive features of atom interfer¬ 
ometry sensors is the ability to perform differential ac¬ 
celeration measurements by simultaneously interrogat¬ 
ing two separated atomic clouds with high rejection of 
common-mode vibration noise, as demonstrated in grav¬ 
ity gradiometry applications (3| 13 • I 11 principle, such a 
scheme can be extended to an arbitrary number of sam¬ 
ples, thus, providing a measurement of higher-order spa¬ 
tial derivatives of the gravity field. Geophysical models of 
the Earth’s interior rely on the inversion of gravity and 
gravity gradient data collected at or above the surface 
[25j. The solution to this problem, which is, in general, 
not unique, leads to images of the subsurface mass distri¬ 
bution over different scale lengths 0]. In this context, 
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the simultaneous determination of gravity acceleration 
and its derivatives improves the inversion procedure by 
introducing additional constraints for the valid solutions. 
Gravity gradient surveys are already used to detect short- 
wavelength density anomalies or in situations where the 
vibration noise seriously limits absolute gravity measure¬ 
ments. The second derivative of the gravity field can vary 
by several orders of magnitude when measured across 
shallow density anomalies, promising high spatial resolu¬ 
tions and sharp signals for their localization [0 . Simul¬ 
taneous in situ measurements of the gravity acceleration 
and its derivatives can also be used for remote sensing 
to estimate the evolution of the gravitational field along 
the direction of the local plumb line. Such a method 
could find interesting applications in regional height sys¬ 
tems to measure differences in the gravitational potential 
with respect to a reference station, e.g., located on the 
geoid 0]. Indeed, in the presence of shallow density 
anomalies, the knowledge of both the gravity gradient 
and the curvature can provide centimeter-level resolu¬ 
tion (~ 0.1 m 2 /s 2 ) in the measurement of differential 
geopotential heights by integrating the gravity field over 
baselines of several hundreds of meters. The simulta¬ 
neous measurement of gravity gradient and higher-order 
derivatives would also help with correcting for Newtonian 
noise in future gravitational-wave detectors [0. 


In this Letter, we report for the first time the direct 
measurement of the gravity-field curvature generated by 
nearby source masses, as suggested in Ref. [9|. Our atom 
interferometer, which simultaneously probes three freely 
falling samples of 87 Rb, is able to perform measurements 
of gravity, gravity gradient, and curvature along the ver¬ 
tical direction at the same time, opening new perspec¬ 
tives for geodesy studies and Earth monitoring applica¬ 
tions. Using this scheme, we also demonstrate a new 
method to measure the Newtonian constant of gravity. 


The details of the experimental apparatus can be found 
in Refs. 00 . In the following, a description of the 
measurement sequence and data analysis will be pro¬ 
vided, with particular emphasis on the new features in- 
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troduced by the third atomic sample and the gravity cur¬ 
vature determination. 

A magneto-optical trap (MOT) with beams oriented 
in a 1-1-1 configuration collects 87 Rb atoms and launches 
them vertically at a temperature of about 4 fjK. A high- 
flux source based on a 2D MOT provides large atom num¬ 
bers (~ 10 9 ) in short loading times (~ 40 ms). Larger 
atom numbers could be obtained by using the juggling 
technique [30]; however, only a direct launch can be read¬ 
ily implemented in our measurement cycle and extended 
to three or more samples. We launch three atomic clouds 
along the vertical direction separated by ~ 30 cm, which 
reach the apogees of their atomic trajectories simulta¬ 
neously. A series of velocity selection and blow-away 
pulses prepares the samples in the magnetically insen¬ 
sitive l-F = 2, m.F = 0) sublevel of the 87 Rb ground 
state. The Mach-Zehnder interferometer simultaneously 
addresses the three clouds with a7r/2 — 7r — 7 t/2 sequence 
of vertical velocity-selective Raman pulses |3l| • The Ra¬ 
man lasers, with effective wave vector k e g- ~ 16 x 10 6 
m , are resonant with the 6.8 GHz two-photon transi¬ 
tion |F = 2 ,rriF = 0) —► |F = 1, To/? = 0) of the 87 Rb 
ground state and have a 2 GHz red detuning with respect 
to the 5 2 S 1/2 |.F = 2) —> 5 2 P 3 / 2 | F = 3) transition to the 
excited state. The sequence has a duration of 2 T = 320 
ms. The 7r pulse lasts 24 /zs and occurs 5 ms after the 
atomic clouds have reached their apogees. The interfer¬ 
ence fringes are obtained by measuring the normalized 
population in one of the two hyperfine levels of the 87 Rb 
ground state. We use a set of high-density source masses, 
for a total of 516 kg, to enhance the gravity-field curva¬ 
ture sensed by the three atomic samples. The source 
masses are composed of 24 tungsten alloy (INERMET 
IT180) cylinders [32]. They are positioned on two tita¬ 
nium platforms and distributed in hexagonal symmetry 
around the axis of the interferometer tube. The verti¬ 
cal position of the platforms is accurately controlled by 
precision screws synchronously driven by stepper motors 
and measured by an optical readout system. The ex¬ 
perimental setup is shown in Fig. [I] together with the 
axial acceleration profile due to the source masses and 
the Earth’s gravity. 

Atomic gravity gradiometers use the same Raman 
lasers to simultaneously probe two spatially separated 
atom clouds on the same interferometric sequence. In 
this configuration, vibration noise that couples into the 
phase of the Raman lasers is seen as common mode 
and can be efficiently rejected. As a consequence, when 
the normalized atomic populations (x, y ) simultaneously 
measured at the two spatially separated interferometers 
are plotted in 2D space, an ellipse is obtained. Common¬ 
mode phase noise affecting the fringes of the two atom in¬ 
terferometers distributes the experimental points around 
the ellipse, whose shape carries information on the grav¬ 
ity gradient between the two clouds [33|. We extend this 
idea by introducing a third atom interferometer. In this 
case, the normalized atomic populations ( x , y, z ) mea¬ 
sured at the output ports of the three atom interferom¬ 



FIG. 1. (color online) a) Scheme of the experiment. s 'Rb 
atoms are trapped and cooled in a MOT. Three atomic clouds 
are launched in rapid sequence along the vertical direction 
with a moving optical molasses. Near the apogees of the 
atomic trajectories, a measurement of the vertical accelera¬ 
tion sensed by the three clouds is performed by Raman inter¬ 
ferometry. External source masses are positioned in order to 
maximize the average gravity curvature at the three clouds’ 
positions, b) Gravitational acceleration along the symme¬ 
try axis (o 2 ) produced by the source masses and the Earth’s 
gravity gradient; a constant value accounting for the Earth’s 
gravitational acceleration was subtracted. The spatial regions 
of the three atom interferometers are indicated by the thick 
red lines. 


eters are distributed around a Lissajous ellipse lying in 
3D space. 

The ellipse best fitting the data points in 3D space can 
be expressed with the parametric equations 


{ x(6) = A + B sin# , 

y{9) = C + D sin(9 + ipi), (1) 

z(9) = E + Fsin(9 + ipi + v? 2 ) • 


Here A, B , C , D , E , and F represent the amplitude and 
offset of the fringes of the three atom interferometers, 9 
is the phase angle parameter, which varies randomly due 
to common-mode vibration noise, and and <p 2 are the 
phase angles proportional to the differential accelerations 
between adjacent interferometers. Fitting an ellipse to 
points in 3D space can be recast as a 2D problem 34]. 
Given a set of n data points ( Xi,yi,Zi ), the y 2 function 
can be written as 


X 2 OC XMeUipse^i, Vi, *i) + *■?] , 

2=1 


( 2 ) 


where ellipse denotes the Euclidean 2D distance be¬ 
tween the ellipse and the projection of the point on the 
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plane of the ellipse, and hi is the 3D point-plane dis¬ 
tance. We assume equal uncertainties on all experimen¬ 
tal points. The y 2 function is then evaluated and min¬ 
imized with respect to the eight parameters of Eq. [I] 
This approach can be easily generalized to N interfer¬ 
ometers (aq, xjv) for the measurement of higher-order 
derivatives of the gravity field along the vertical axis. 

It is worth pointing out that adding an extra dimension 
(i.e., a third atom interferometer) opens the possibility 
of accurately measuring small gradiometric phase shifts 
introducing negligible bias on the fit results. In a two- 
cloud configuration and in the presence of a small gravity 
gradient, the two output fringes are almost in phase and 
the ellipse degenerates to a line. On the other hand, in 
the presence of a third cloud, even if ip 2 ~ 0 , ipi can 
be made quite large, e.g., by pulsing a magnetic field 
at the location of the third interferometer. In this case, 
|y >2 — fi\ ^*0 and, even in presence of noisy data, p >2 can 
be reliably extracted from the fit in 3D space. Figure 
[ 2 ] compares the bias errors introduced by three different 
methods: a least-square fit in 2D, a least-square fit in 3D 
space, and a Bayesian analysis, all as a function of ip 2 , 
when ipi is kept at n/2. Simulated data are affected by 
significant Gaussian noise at detection (ad = 0.01) and 
present a fringe contrast of 0.3. The plot clearly shows 
how the third interferometer becomes instrumental for 
precision gravity gradiometry. The bias errors introduced 
by the 2D fit are non-negligible. In addition, for (f >2 < 
0.1, the 2D fit routine fails to converge. The Bayesian 
analysis performs better than the 2D fit, but it introduces 
significant biases at small phase angles when the a priori 
knowledge of the noise affecting the data varies by 10 %. 
On the other hand, the 3D fit is very robust and, in 
contrast with the Bayesian method, does not require any 
a priori knowledge of the noise on the experimental data. 


Our setup has been used to perform a direct mea¬ 
surement of the gravity-field curvature generated by the 
source masses. One of the most critical aspects of the 
measurement is the presence of spurious and nonhomo- 
geneous magnetic fields in the interferometer region. Be¬ 
cause of the spatial separation between the three atomic 
clouds, primarily imposed by the MOT loading time, the 
gravity curvature measurement is averaged over a total 
distance of about 60 cm. In this configuration, the lower 
and the upper atomic samples are close to the edges of the 
p -metal shield surrounding the vertical tube, where the 
passive attenuation of external magnetic fields is lower 
and the internal bias held is less homogeneous. To re¬ 
duce this source of systematic errors, the sign of the 
effective wave vector k e g- of the Raman lasers is peri¬ 
odically reversed during data acquisition 35] - This is 


achieved by selecting a different Raman counterpropa- 
gating beam pair by properly adjusting the frequency 
detuning to compensate for the Doppler shift induced by 
the atomic motion in the gravity field. In this way, phase 
shifts that do not depend upon the effective wave vec¬ 
tor, e.g., second-order Zeeman shifts or ac Stark shifts, 
are rejected when taking the difference between measure- 
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FIG. 2. (color online) Bias error in the differential phase 
from the 2D elliptical fit (red circles), the new 3D fit rou¬ 
tine (black squares), and the Bayesian method for different 
<P 2 angle values; ipi is kept fixed at tt/2. Synthetic data are 
generated with a Gaussian detection noise (ad = 0.01) and a 
fringes contrast of 0.3. In the Bayesian analysis, we feed the 
algorithm both with the exact Gaussian noise as used in the 
simulation (blue triangles) and with the value obtained after 
introducing a 10% error on ad (green rhombi). For <^2 <0.1 
rad, the numerical minimization fails in the 2D fit. The 
knowledge of the noise model affecting the data becomes crit¬ 
ical in the Bayesian analysis when the phase angle approaches 
zero. 


ments performed with opposite fc e ff- Submillimetric verti¬ 
cal overlap of the interferometer’s arms has been achieved 
by properly adjusting the Raman frequency ramp and the 
timing of the velocity selection. Their transverse overlap 
is ensured by using the same launch sequence. The hori¬ 
zontal velocity spread due to the finite transverse atomic 
temperature is expected to introduce noise and system¬ 
atic shifts on the ellipse phase angle via the Coriolis accel¬ 
eration. Because of the double differential nature of the 
gravity-field curvature measurement, the effect of Cori¬ 
olis accelerations depends on the difference between the 
relative initial velocities of the atomic clouds at the two 
adjacent gravity gradiometers. To further reduce this ef¬ 
fect, the mirror retroreflecting the Raman laser beams is 
rotated to compensate for the Earth’s rotation 36. ; 37| . 
Two data sets of 720 points (2.5 s of measurement time 
per point), one for each of the two fc e ff opposite directions 
(t and j. ), have been collected and analyzed. Figure [3] 
shows a typical plot of the data points measured at the 
three conjugated atom interferometers, together with the 
ellipse in 3D space best fitting the data. The values for 
ipi and ip 2 are given by 

V?i = (<m “ m )/ 2 > V 2 = (v? 2 ,f - <P 2 , 4)/2 . (3) 

From the measurement of the clouds’ separation, d = 
(0.3098 ± 0.0002) m, it is possible to evaluate the av¬ 
erage gravity gradients 71,2 = ^ 1 , 2 / (dfc e ffT 2 ), obtaining 
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FIG. 3. (color online) Typical three-dimensional Lissajous 
figure obtained by plotting the output signal of the upper 
atom interferometer as a function of the lower and central 
one (red circles) and ellipse in 3D best fitting the data (black 
line). Orthogonal projections on the three Cartesian planes 
are also shown. 


71 = (-4.112 ± 0.008) x 10 ~ 6 s -2 and 72 = (0.223 ± 
0.003) x 10 ” 6 s -2 , and, thus, the average gravity curva¬ 
ture (f = (1.399 ± 0.003) x 10 -5 s _ 2 m _1 . This measure¬ 
ment is consistent with the value £ s i m = 1-397 x 10 -5 
s _ 2 m _1 obtained from our Monte Carlo model [38| . 
which accounts for the source masses and additional con¬ 
tributions in the immediate vicinity of the atomic clouds. 

The measurement of the second derivative of the grav¬ 
ity acceleration is also an interesting tool for determining 
the Newtonian gravitational constant G, as proposed in 
Ref. [1^. The method consists of performing two si¬ 
multaneous gravity gradient measurements in the pres¬ 
ence of heavy source masses. The Earth’s gravity gradi¬ 
ent contribution is rejected when calculating the differ¬ 
ence between the two measurements without any need 
for modulating the position of the masses. In this way, 
systematic effects introduced by deformations and tilts 
of the structure holding the masses can be removed. For 


such an experiment, it becomes important to optimize 
the distribution of the source masses to generate three 
quasistationary regions to host the conjugated atom in¬ 
terferometers, thus, reducing the systematics arising from 
the positioning errors of the atomic clouds. Even if not 
specifically designed for this purpose, we used our appa¬ 
ratus to perforin a proof-of-principle experiment. With 
the source masses and the atomic clouds positioned as in 
Fig. 1, we measured ^Veas = ip 2 — ¥7 and compared it 
with cE» s i m obtained from our single-particle Monte Carlo 
simulation. We obtained 4> meas = (0.5533 ± 0.0006) rad, 
which is in good agreement with <l> s i m = 0.5528 rad. The 
short-term sensitivity of 3.8 x 10“ 1 2 3 G at 1 s is comparable 
with the one obtained in Refs. [ 1 , EH, Hl| by alternating 
the source masses position. An extensive evaluation of 
the systematic error sources that are affecting the mea¬ 
surement is beyond the scope of this work. 

In conclusion, by using three simultaneous atom in¬ 
terferometers, we have measured for the first time the 
component of the gravity curvature produced by nearby 
source masses along one axis. The new analysis method 
based on an elliptical fit in 3D space has proven to be 
very robust with respect to amplitude noise and immune 
from noise-induced systematic shifts. The scheme has 
also been used to perform a proof-of-principle measure¬ 
ment of the Newtonian gravitational constant based on 
two simultaneous gravity gradient measurements. Sensi¬ 
tivity and long-term stability of the G measurement are 
comparable with our previous work, opening the possi¬ 
bility for further improvements after optimization of the 
distribution of the masses and the position of the atomic 
clouds. This method can be extended to multiple inter¬ 
ferometers with small spatial separation (~ 5 — 10 cm) in 
order to reconstruct acceleration profiles with high reso¬ 
lution and measure higher-order derivatives of the grav¬ 
itational acceleration. 
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